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INTRODUCTION 


1 .  ORIGIN  OF  PROGRAM 

This  program  was  developed  in  The  Hydrologic  Engineering  Center  by 
R.  G.  Willey  and  H.  E.  Kubik.  A  two-part  FORTRAN  II  version  was  written 
in  1968;  a  September  1970  version  was  written  in  FORTRAN  IV  and  included 
increased  capability,  and  this  version  basically  modifies  the  input  to 
conform  with  present  standards. 

It  is  recognized  that  nearly  all  computer  centers  have  access  to  multiple 
linear  regression  programs,  but  this  program  includes  some  different 
statistical  philosophy  and  some  capability  that  is  particularly  useful  for 
regional  analyses. 

2.  CAPABILITIES  OF  PROGRAM 

The  program  makes  a  multiple  linear  regression  analysis  and  has  the 
following  special  features: 

a.  Automatic  deletion.  The  analysis  is  first  made  with  all  of  the 
specified  independent  variables  3nd  then  the  least  significant  variable  is 
deleted  and  the  analysis  repeated. 

b.  Variable  selection.  All  independent  and  dependent  variables  can 
p-  stored  and  then  only  those  variables  desired  for  a  particular  analysis 
"'ay  be  used. 

c.  Combination  of  variables.  New  variables  may  be  computed  from 

the  input  variables,  e.g.,  a  parameter  A/5  [basin  area (A)  times  the  square 
root  of  main  channel  slope  (S)]  can  be  computed  when  area  and  slope  have 
been  input. 

d.  Transformations .  The  variables  may  be  transformed  (square  root, 
logarithmic  or  reciprocal)  to  more  nearly  linearize  the  relations. 


‘r 

ft 


limi1  8  j 


e.  Input  regression  parameters.  The  regression  constant  and 
regression  coefficients  may  be  input  and  the  residuals  (difference 
between  observed  and  calculated),  the  adjusted  multiple  determination 
coefficient  and  the  standard  error  of  estimate  will  be  computed. 

3.  HARDWARE  AND  SOFTWARE  REQUIREMENTS 

This  program  has  been  developed  and  tested  primarily  on  the  UNIVAC  1108 
and  the  CDC  7600  computer.  No  tape  units  are  required. 


PROGRAM  DESCRIPTION 


4.  PROGRAM  ORGANIZATION 

The  program  consists  of  a  main  program  and  two  subroutines.  Subroutine 
COMB  is  used  to  compute  a  new  variable  based  on  a  combination  of  one  or 
more  variables.  Subroutine  CROUT  solves  a  set  of  simultaneous  equations 
to  obtain  the  regression  coefficients.  The  proqram  will  exit  normally 
only  in  the  main  program.  A  macro-flow  chart  of  the  program  organization 
is  shown  in  figure  1. 

5.  THEORETICAL  ASSUMPTION  AND  LIMITATIONS 

The  computations  are  made  in  accordance  with  the  procedures  given  in 
"Statistical  Methods  in  Hydrology"  by  Leo  R.  Beard,  January  1962.  The 
independent  variables  are  deleted,  in  turn,  based  on  the  minimum  adjusted 
partial  determination  coefficient  (r  ^). 

The  adjusted  multiple  determination  coefficient  is  used  to  compute  the 
standard  error  of  estimate.  Care  must  be  taken  in  insuring  that  there 
are  more  observations  than  variables.  If  the  number  of  observations  are 
very  close  to  the  number  of  variables,  the  results  may  be  unreasonable. 

Variables  with  zero  or  very  little  (.000001)  variance  are  deleted  from 
the  analysis.  A  diagnostic  is  printed  stating  the  name  of  the  variable 
deleted  because  of  this  limitation. 

6.  METHODS  OF  COMPUTATION 

The  regression  analysis  is  performed  on  either  the  transformed  or  non- 
transformed  variables  in  a  like  manner.  The  basic  equation  is: 


where  the  X-j  values  are  the  independent  variables  and  the  Y  value  is 
the  dependent  variable,  "n"  is  the  number  of  independent  variables  in 
the  analysis.  The  regression  coefficients  (b-j)  are  calculated  using  the 
Crout  Method  (Exhibit  1)  for  solution  of  the  following  equations: 

b-,  A  :  (x]x2)  b2  +  •  •  ■  +  '  (x.1xiJ  t»n  =  "(yx-j ) 
t ( X-j x.> )  b,  +  t(x2)“  b2  +  •  •  •  T(x2xn)  tr  --  t(yx2) 


(2) 

(xlxn)  bl  +  '(YnJ  b2  +  •  *  •  +  :(xn)2  bn  =  :(yxn} 


where  x.  is  the  deviation  of  the  value  X  for  variable  i,  the  observation 
numb01'-  ^°ing  implied,  from  the  mean  (X.).  T.(x.  )2  and  f(x^x.)  can  be 

determined  from  the  following  equations:  J 

::(xi)2  =  (Xi  -  X.)2  =  ::(Xi  )2  -  (rx.  )2/N  (3) 

:-(x.Xj)  -  :-(X.  -  X.)(Xj  -  X.)  =  ^(X.Xj)  -  "X.  T. X j / N  (4) 


where  N  is  the  number  of  events  or  observations. 

The  regression  constant  (a)  is  calculated  by  use  of  the  equation: 

a  =  Y  -  X  b.X.  (5) 

where  Y  is  the  mean  of  the  dependent  variable  for  N  events  and  the  X. 
values  are  the  means  of  each  of  the  independent  variables  for  N  events. 

2  ~2 

The  unadjusted  (R  )  and  the  adjusted  (R  )  determination  coefficients 
are  determined  by  the  following  equations: 

P  b-  2(yx,)  +  bp  r. ( yxp )  +  .  .  .  +  bn  r.(yxn) 

R  =  - - - - - -  y-  - - - — -  (6) 

-(yr 
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=  1  -  (1  -  R2)(N  -  1 )/ (N  -  NVAR) 


(7) 


where  NVAR  is  total  number  of  variables  (both  dependent  and  independent) 
in  the  analysis. 


The  partial  determination  coefficient  for  a  variable  is  the  percent 
increase  in  unexplained  variance  caused  by  deleting  that  variable  from 
the  regression  equation.  The  adjusted  partial  determination  coefficient 
is  calculated  by  the  following  equation: 


—2  (1-R  •]  23)  ~  0"^  ]  2 34 ^ 

r  14 .23  =  - - — - Lz£2L 

1  -  R\ 

1  .  2j 


where  the  first  subscript  for  the  adjusted  partial  determination 
coefficient  indicates  the  dependent  variable,  the  second  indicates  the 
independent  variable  whose  partial  determination  coefficient  is  being  com¬ 
puted,  and  tiie  subscripts  after  the  decimal  indicate  the  other  independent 
variables  involved  in  the  computation. 

The  standard  error  of  estimate  (S  )  is  calculated  by  the  following  equation 


3e  =  >'  (1-R2)':(y)2/(N-1)  =  Sy  /l-R2 


where  S  is  the  standard  deviation  of  the  dependent  variable. 

y 


When  the  residual  prediction  option  is  specified  the  "RESIDUAL"  column 
is  calculated  as  the  difference  Ketween  the  observed  and  predicted  values. 
The  "ERROR  AS  DEVIATE"  column  is  calculated  as  the  residual  divided  by 
the  standard  deviation.  The  "RATIO"  column  is  calculated  as  the  residual 
divided  by  the  observed  value. 
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EXHIBIT  1 
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One  of  the  best  methods  for  solving  systems  of  linear  equations 
on  desk  calculating  machines  was  developed  by  P.  D.  Crout  in  1941. 

This  method  is  based  on  the  elimination  method,  with  the  calculations 
arranged  in  systematic  order  so  as  to  facilitate  their  accomplishment 
on  a  desk  calculator.  In  this  method  the  coefficients  and  constant 
terms  of  the  equations  are  written  in  the  form  of  a  "matrix,"  which  is 
a  rectangular  array  of  quantities  arranged  in  rows  and  columns. 


The  method  is  beat  explained  by  an  example.  Suppose  that  in  a 
multiple  correlation  analysis  it  is  required  to  solve  the  following 
system  of  linear  equations  to  obtain  the  unknown  values  of  b  ,  b,,  b. 
and  b  .  5 
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A  continuous  oheok  on  the  computations  as  they  progress  may  be  obtained 
by  adding  to  the  matrix  of  the  abovs  system  a  column  of  u's,  such  that 
u-p+q+r+e+t.  The  matrix  and  oheok  oolumn  are  written  as 
follows* 


The  Crout  Method  was  presented  at  the  AIEE  Summer  Convention  in  June  1941 
by  Prescott  D.  Crout.  The  method  was  developed  by  Gauss  and  refined  by 
Duo  1 i 1 1 1 e . 
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The  elements  p^,  q^,  r^  and  3^  form  the  "principal  diagonal"  of 

the  matrix.  Examination  of  the  original  equitions  shows  that  the  co¬ 
efficients  are  symmetrical  about  the  principal  diagonal,  i.e.,  «  p0 , 

rl  *  ?3’  r2  "  q3’  S1  “  ps-’  s2  “  c-4*  “na  s5  “  r4* 

This  is  characteristic  of  the  system  of  equations  to  be  solved  in  any 
multiple  correlation  analysis.  Because  of  this  symmetry,  the  computa¬ 
tions  are  considerably  simplified.  While  the  Crout  method  may  be  used 
to  solve  any  system  of  linear  equations,  the  computational  steps  given 
here  are  applicable  only  to  those  with  symmetrical  coefficients. 

The  solution  consists  of  two  parts,  viz.,  the  computation  of  a 
"derived  matrix"  and  the  "back  solution."  Let  the  derived  matrix  be 
denoted  as  follows: 
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The  elements  of  the  derived  matrix  are  computed  as  follows* 


?1  ■  P1  P2  *  p2  ?3  •  p5  P4  p4 

h  ■  „  r,  s ,  t.  u 

■i  -  —  R1  -  —  S  -  _i  T  -  JL  u  -  -1 

■4  Pt  n  J4r>  A  n.  A  D. 


‘  42  '  ?241  '}  '  q;  ?34i  R2  " 


X  '  “4  '  P4Q1  S2  ■  P2 


j.  n  T) 

v2‘nP2 


R3  “  r5  ”  SR2  "  P3R1  \  "  r4  "  V*2  "  PdRi 
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s,  -  — 

5  R3 


_  t3.y3.T1P, 

3  *3 


B_  -  “3-V3-V3 
5  S3 


S4  -  84  -  R^Sj  -  ft4S2  -  P4Sj 

VT3VT2VT1P4  _  U4“D3R4~R2^4-R1P4 

m  . - «: — -  0.  -  —  <• — -■ -  -  — 


The  general  pattern  of  the  above  computations,  which  may  be  applied 
to  a  system  containing  any  number  of  equations,  is  as  follows t 

(1)  The  first  column  of  the  derived  matrix  is  copied  from  the  first 
column  of  the  given  matrix. 

(2)  The  remaining  elements  in  the  first  row  of  the  derived  matrix 
are  computed  by  dividing  the  corresponding  elements  in  the  first  row  of 
the  given  matrix  by  the  first  element  in  that  row. 

(3)  After  completing  the  n—  row,  the  remaining  elements  in  the  (n+l) — 
column  are  computed.  Suoh  an  element  (X)  equals  the  corresponding  element 
of  the  given  matrix  minus  the  product  of  the  eloment  immediately  to  the  loft 
of  (X)  by  the  element  immediately  above  the  principal  diagonal  in  the  s .0 
column  as  (X),  minus  the  product  of  the  second  element  to  the  left  of  (X) 

by  the  second  element  above  the  principal  diagonal  in  the  same  column  Qr.  (X), 
eto.  After  each  element  below  the  principal  diagonal  is  recorded,  and  while 
that  element  is  still  in  the  calculator,  it  is  divided  by  the  element  of 
the  principal  diagonal  which  is  in  the  same  column.  The  quotient  is  the 
element  whose  location  is  symmetrical  to  (X)  with  respeot  to  the  principal 
diagonal, 
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(4)  When  the  elements  in  the  (n+l)-^  column  and  their  symmetrical 
counterparts  have  been  recorded,  the  (n-rl)iil  row  will  be  complete  except 
for  the  last  two  elements,  which  are  next  computed.  Such  an  element  (X) 
equals  the  corresponding  element  of  the  given  matrix  minus  the  produot 
of  the  element  immediately  above  (X)  by  the  element  immediately  to  the 
left  of  the  principal  diagonal  in  the  same  row  as  (X),  minus  the  product 
of  the  second  element  above  (X)  by  the  second  element  to  the  left  of  the 
principal  diagonal  in  the  same  row  as  (X),  eto.,  all  divided  by  the 
element  of  the  principal  diagonal  in  the  same  row  as  (X). 


The  check  column  (U)  of  the  derived  matrix  serves  as  a  continuous 
check  on  the  computations  in  that  each  element  in  the  oolumn  equals 
one  plus  the  Siam  of  the  elements  in  the  same  row  to  the  right  of  the 
principal  diagonal.  That  is, 


1  + 


S  +  T 
1  1 


U„ 


i  +  r2  +  s2 


-t-  r. 


1  ♦  S5  +  Tj 


This  check  should  be  made  after  completing  each  row. 

The  elements  of  the  derived  matrix  to  the  right  of  the  principal 
diagonal  form  a  system  of  equations  which  may  now  be  used  to  compute 
the  unknown  values  of  b2>  b^,  and  b^  by  successive  substitution. 

This  is  known  as  the  "back  solution."  The  computations  are  as  follows* 
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It  is  very  important  that  the  computations  be  carried  to  a  suffi¬ 
cient  number  of  digits,  both  in  computing  the  coefficients  and  constant 
terms  of  the  original  equations,  and  in  computing  the  elements  of  the 
derived  matrix.  It  is  possible  for  relatively  small  errors  in  the 
coefficients  and  constant  terms  of  the  original  equations  to  result  in 
relatively  large  errors  in  the  computed  solutions  of  the  unknowns.  The 
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greatest  source  of  error  in  computing  the  elements  of  the  derived  matrix 
arises  from  the  loss  of  leading  significant  digits  by  subtraction.  This 
must  be  guarded  against  and  can  be  done  bv  carrying  the  computations  to 
more  figures  than  the  data.  A3  a  general  rule,  it  is  recommended  that 
the  coefficients  and  constant  terms  of  the  original  equations  be  carried 
to  a  sufficient  number  of  decimals  to  produce  at  least  five  significant 
digits  in  the  smallest  quantity,  and  that  the  elements  of  the  derived 
matrix  be  carried  to  one  more  decimal  than  thi3,  but  to  not  less  than 
six  significant  digits. 
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.  75  5 

.  168 

.413 

1 . 039 

DT 

19  4  7 

.  960 

.  280 

.410 

.979 

ED 

NT 

1 

77 

DATA 

FROM  PREVIOUS 

ANALYSIS 

TR 

7 

0 

0  1  1 

NT 

71 

res 

T  NO.  2 

JANUARY  1975 

73 

DAT 

A  WITH 

SPECIAL 

FORMAT 

T7 

COMPUTATION 

OF  NEW 

VARIABLE 

J1 

4 

1 

1 

NM 

AREA  LENGTH 

SLOPE 

LOG  Q 

F7 

(A3,  A- 

5,  7FB.0/9X, 

F8. 0) 

C  T 

9 

1  L/S*  *  .  5 

CO 

-3 

1 

;3 

5 

6 

cc 

.  5 

rpp 

4 

7 

0  0  17 

DT 

TDP  R 

230 

38.8 

1  .  2? 

DT 

I  DP  R 

3.082 

DT 

BUF  C 

19.4 

10.9 

15.42 

or 

BUF  C 

2.  377 

or 

MCD  C 

7.52 

7.0 

9.66 

DT 

MCD  C 

2.062 

DTSALT  1 

32.5 

10.9 

13.04 

DTSALT  1 

2. 568 

DTS ALT  2 

1  14 

CT 

vT 

CO 

3.01 

DTSALT  2 

3. 004 

DT 

7  DP  R 

63  5 

98.  1 

1.06 

DT 

2DP  R 

3.557 

ED 

ED 

MJ 

T 1  TEST  NO.  3  JANUARY  1975 

T2  USE  DATA  FROM  TEST  NO.  1 

T3  READ  REGRESSION  COEFFICIENTS 


J1 

4 

i 

NM 

LOG  Q 

LOG  SNO 

LOG  GWLOG 

PRCP 

rp 

1  ' 

1111 

RP 

-.22 

1 . 6 

1 . 0 

.  3 

DT 

1936 

.939 

.  399 

.  325 

.710 

DT 

1937 

.  94  5 

.  343 

.  385 

.634 

DT 

1938 

1. 052 

.  369 

.408 

.  886 

DT 

1979 

.744 

.  246 

.  428 

.531 

DT 

1 94  0 

.666 

.  181 

.316 

1  .  027 

DT 

1941 

1 . 091 

.  297 

.460 

1.315 

DT 

194"’ 

1 . 060 

.  299 

.511 

1 . 097 

or 

194  3 

.  892 

.  354 

.  379 

.70  7 

DT 

1944 

1.021 

.  295 

.  395 

1  .  240 

Dr 

194  5 

.  920 

.  321 

.  7  76 

1.091 

DT 

194  6 

.  75  5 

.  168 

.413 

1 . 038 

ITT 

194J 

.  960 

.  ">80 

.410 

.979 

ED 


test  data  ultitt 


******************************** 

MULTIPLE  LINEAR  REGRESSION 

704-G1-L20J0  JAN  1975 

REVISED  7  9N  1933 

******************************** 


TEST  NO.  1  JANUARY  1^75 
STATISTICAL  METHODS  { BEARD  1962)  EX  33 
FIRST  ANALYSIS-DATA  SUPPLIED 


NVAL 

NO  1MB 

TPRNT 

TFRMT 

!  DELE 

4 

0 

0 

0 

0 

DELTA 

0.0000 


+  +  4-  t-  ANALYSIS  NO  1 
DEPENDENT  VARIABLE  -- 


IDRP 

.  .  .  .  .TRANSFORMATION 

CODES . 

NOBR 

1 

1  ' 

1  1  1  0  0  0  0  0  0  0 

0  0  0  0  0  0  0 

0 

VARIABLE 

TRANSFORMATION 

LOG  Q 

NONE 

LOG  3 NO 

NONE 

LOG  GW 

NONE 

LOG  PRCP 

NONE 

INPUT  DATA 

OBS  NO  OB 

S  ID 

LOG  O 

LOG  SNu 

LOG  GW 

LOG  PRCP 

1 

1 936 

0.939 

0  .  399 

0.325 

0.710 

2 

1937 

0 . 94  3 

0.343 

0.335 

0.634 

3 

1938 

1 .032 

0.369 

i).  403 

0 . 836 

4 

193  9 

0.744 

0 . 246 

0.423 

0.531 

3 

194  0 

3 .666 

0.131 

0.316 

1 .027 

A 

1  94  1 

1.031 

0 . 297 

0 . 460 

1.315 

7 

194  2 

1 . 060 

0 . 299 

0.511 

1 . 097 

3 

1 94  3 

7 . 302 

0. 354 

0.779 

0.707 

? 

194  4 

1.0O1 

0 . 195 

0.  395 

1 .  240 

10 

194  3 

0.920 

0 . 321 

0.376 

1.091 

1  1 

1  34  6 

0.75  3 

0.  168 

0.413 

1  .  038 

1  2 

1947 

0 . 960 

0. 230 

0.410 

0.979 

I  Mi 


2  o  I  1 


4-4-4- 

LOG  n 

IFORC 

0 


o'i'A T  !  ST  1  OS  OF 

DATA 

BASED  ON  12 

OBSERVATIONS 

STANDARD 

VARIABLE 

AVERAGE 

VARIANCE 

DEVIATION 

LOG  SNO 

0.2360 

0 . 0050 

0.0704 

LOG  GW 

0 .4005 

0.0023 

0.0531 

;PG  PROP  0.9421  0.0572  0.2302 

Q  0.0105  0.0181  0.1346  DEPENDENT  VARIABLE 


SIMPLE  CORRELATION  COEFFICIENTS 


VART  ABLE 

LOG  SNO 

LOG  GW 

LOG  PRCP 

LOG  Q 

LOG  SNO 

1  . 0000 

0.0000 

-0.0459 

0.6309 

LOG  GW 

0.0000 

1 .0000 

0.1275 

0.4170 

LOG  PRCP 

-0.0459 

0.1275 

1.0000 

0.2011 

LOG  Q 

0.6303 

0.4170 

0.2011 

1.0000 

INDEPENDENT  REGRESSION  PARTIAL 

VARIABLE  COEFFICIENT  DETERMINATION 

COEFFICIENT 

LOG  SNO  1.521805  0.9105 

LOG  GW  1.01 2931  0.6314 

LOG  PRCP  0.273387  0.7451 


STANDARD 

REGRESSION  R  R  BAR  ERROR  OF 

CONSTANT  SQUARE  SQUARE  ESTIMATE 

-0.223704  0.9437  0.9226  0.0374 
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VARIABLE  DELETED  IS 


LOG  GW 


INDEPENDENT 

REGRESS  ION 

PARTIAL 

VARIABLE 

COEFFICIENT 

DETERMINATION 

COEFFICIENT 

LOG  SNO 

1.64731 3 

0.7467 

LOG  PROP 

0 . 3437  3  0 

0. 5063 

REGRESSION 

R 

R  BAR 

STANDARD 

ERROR  OF 

CONSTANT 

SQUARE 

SQUARE 

ESTIMATE 

0.103475 

0.3011 

0 . 7570 

0.0663 

DEPENDENT 

VARIABLE  LOG 

o 

ft. 

DBS  NO  OBS  ID 

OBSERVED 

CALCULATED 

RESIDUAL 

ERROR  AS 

DEVIATE 

RATIO 

1 

1  3  3  5 

3.93  9 

1  .  303 

-0.069 

-0.515 

-0 . 074 

2 

10  3  7 

0 . 94  5 

0.390 

0.05  5 

0.41  .3 

0.0  59 

T 

1338 

1.05  3 

1 . 020 

0.033 

0 . 2  36 

0 . 030 

4 

133 '3 

0.744 

0.711 

0.033 

0.  243 

0.044 

5 

1 94  0 

0.666 

0 . 760 

-0.094 

-0.697 

-0. 141 

6 

1941 

1.031 

1 . 051 

0.030 

0.221 

0.027 

7 

1 94 .7 

1 . 060 

0.979 

0.031 

0.60  5 

0.077 

8 

194  3 

0 . 393 

0.93  3 

-0.041 

-0.306 

-0.046 

9 

194  4 

1.021 

1 . 02  2 

-0.001 

-0.006 

-0.001 

10 

194  5 

0.920 

1.013 

-0.093 

-0.689 

-0.101 

1  1 

1 94  6 

0.755 

0.742 

0.013 

0.095 

0.017 

12 

194  7 

0 . 960 

0.906 

0. 054 

0.400 

0.056 

R  BAR  SQUARED 
STANDARD  ERROR  OF  ESTIMATE 
MEAN  ERROR 
MEAN  SQUARED  ERROR 


0.7570 

0.0563 

0.0000 

0.0033 


t-  4-  V  M 


•  M-  I-  t-  t-  1-4-  I 


4-  t  I-  f  4- 


I’a si1  4  ol  11 


VARIABLE  DELETED  IS  LOG  PRCP 


INDEPENDENT  REGRESSION  PARTIAL 

VARIABLE  COEFFICIENT  DETERMINATION 

COEFFICIENT 


LOG  3 NO 

1 . 2362 1 2 

0. 3979 

STANDARD 

REGRESSION 

R 

R  BAR 

ERROR  OF 

CONSTANT 

SQUARE 

SQUARE 

ESTIMATE 

0 . 533965 

0.4526 

0.3979 

0 .  1044 

ENDS 

NT  VARIABLE  LOG 

Q 

NO 

OBS  ID 

OBSERVED 

CALCULATED 

RESIDUAL 

ERROR  AS 

DEVIATE 

RATIO 

1 

1935 

0.939 

1.052 

-0.113 

-0.340 

-0.120 

2 

193  7 

0.94  5 

0.980 

-0.035 

-0.260 

-0.037 

3 

1938 

1 . 052 

1.013 

0.039 

0 . 286 

0.037 

i 

1939 

0 . 744 

0.355 

-0.111 

-0.827 

-0.150 

5 

194  0 

0.666 

0.772 

-0.  106 

-0.785 

-0. 159 

6 

1941 

1.081 

0.921 

0.  160 

1.  190 

0.  148 

7 

1942 

1 . 060 

0.92  3 

0.137 

1.015 

0.129 

8 

1943 

0.892 

0.994 

-0. 102 

-0.759 

-0.115 

9 

1944 

1.021 

0.918 

0.  103 

0.763 

0.101 

10 

1946 

0.920 

0.952 

-0.032 

-0.236 

-0.034 

1  1 

1946 

0.755 

0 . 756 

0.000 

0.000 

0.000 

12 

1947 

0.960 

0.399 

0.061 

0.453 

0.064 

R  BAR  SQUARED  0.3979 

STANDARD  ERROR  OF  ESTIMATE  0.1044 

MEAN  ERROR  0.0000 

MEAN  SQUARED  ERROR  0.0091 
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*■***-*■**★****  ******************** 

MULTIPLE  LINEAR  REGRESSION 

70  4 -3 1-L202  0  J\N  1075 

REVISE^  JA  I  1933 

***********************  ******** 


TEST  NO.  1  JANUARY  1975 
STATISTICAL  METHODS  ( BEARD  1962)  EX  33 
DATA  FROM  PREVIOUS  ANALYSIS 


+  +  +  t-  ANALYSIS  NO  2  +  +  v  t- 
DEPENDENT  VARIABLE  —  LOG  GW 


I9EP 


. TRANSFORMATION  CODES .  NOBR 

00  1  10  000000000000  0  0 


IRES  IEORC 

0  0 


VARIABLE 

log  2 

LOG  3 NO 
LOG  GW 
LOG  PRCP 


T RANEE  IRMA tion 
NOT  USED 
NOT  USED 
NONE 
NONE 


INPUT  DATA 


NO 

CBS  ID 

LOG  GW 

LOG  PRCP 

1 

1  136 

0.325 

0.710 

2 

1937 

0 . 395 

0.634 

3 

193  8 

0  .  108 

0.896 

4 

1939 

0.428 

0.591 

5 

194  0 

0.316 

1 . 027 

6 

194  1 

0 . 460 

1.315 

7 

194  2 

0.511 

1 . 097 

9 

194  3 

0.379 

0.707 

9 

194  4 

0.395 

1  .  240 

10 

194  5 

0.376 

1 . 091 

1  1 

1946 

0.  113 

1  .  038 

1  2 

194  7 

0 .410 

0.979 

I’.-n-. 


!)  ol  11 


STATISTICS  OP  DATA 


BASED  ON  17 

OBSERVATIONS 

STANDARD 

VARIABLE 

AVERAGE 

VARIANCE 

DEVIATION 

LOG  PRCP 

0.0421 

0.0572 

0. 2392 

LOG  GW 

0.4005 

0. 0020 

0.0531 

DEPENDENT  VARIABLE 

SIMPLE  CORRELATION  COEFFICIENTS 


VARIABLE  LOG  PRCP  LOG  GW 

LOG  PRCP  1 . 1000  0. 1075 


LOG  GW  0.1275  1.0000 


INDEPENDENT 

VARIABLE 

LOG  PRCP 


REGRESSION 

COEFFICIENT 


0.072122 


PARTIAL 

DETERMINATION 

COEFFICIENT 

0.0163 


REGRESSION 
CONSTANT 
0. 332546 


R 

SQUARE 
0. 1057 


STANDARD 
R  BAR  ERROR  OF 

SQUARE  ESTIMATE 

0.0163  0.0526 
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******************************** 

MULTIPLE  LINEAR  REGRESSION 

704-G 1-L2020  JAN  1975 

REVISED  JAN  1993 

******************************** 


TEST  NO. 2  JANUARY  1975 

DATA  WITH  SPECIAL  FORMAT 
COMPUTATION  OF  NEW  VARIABLE 


NVAL  NCOMB  IPRNT  IFRMT  I DELE 

4  10  10 


DELTA 

0.0000 


I 


DATA  FORMAT  (  A2 ,  AS  ,  3F8 . 0  /9X  ,  F9 . 0  ) 


COMBINATION 

1  STEP 

1 

2  3 

4 

5 

6 

OPERATION  CODE 

-3 

1  8 

5 

6 

-2 

CONSTANTS 

1/ 

0.5000 

7  8  9 

1  7  4 


+  +  +  f  ANALYSIS  NO  1  +  *  +•  * 
DEPENDENT  VARIABLE  —  LOG  Q 


IDEP 

4 


. TRANSFORMATION  CODES .  NOBR 

3  0  0  130000000000000  0 


IRES  IFORC 
0  0 


VARIABLE 

AREA 

LENGTH 


TRANSFORMATION 
COMMON  LOG 
NOT  USED 


SLOPE 
LOG  Q 
L/S*  *  .  5 


NOT  USED 
NONE 
COMMON  LOG 


INPUT  DATA 


NO 

OBS 

ID 

AREA 

LOG  Q 

L/S**. 5 

1 

IDP 

R 

230.000 

3.082 

34.429 

2 

BUF 

c 

19.400 

2.377 

2.776 

3 

MCD 

c 

7.520 

2.052 

2.252 

i 

4 

SALT 

1 

32.500 

2 . 558 

3.018 

5 

SALT 

2 

1 14.000 

3 . 004 

21.211 

6 

2DP 

R 

535.000 

3.557 

85.570 

1 
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STATISTICS  OF  DATA 

BASED  ON  A  orSKR'.’AT  t  INS 

ARi  ANTE 

u . jIuj 
0 . 4 526 
0.2936 


>:  a.?:, 


;  *  * 


AS*  *  .  5 

1,00  Q 


A  .  A  rA  ,  h. 


i  .  3 

1.0113 

'>.7750 


STAPLE  CORRELATION  COEFFICIENTS 


VART ABLE 
ARF.A 
L/s**  .  5 
LOG  Q 


AREA 
1 . 0000 
0.90  S3 
0.9936 


L/S** .  5 
0 . 96 S3 
1 . 0000 
0.9566 


LOG  Q 
0.9936 
0.95 66 
1 . 0000 


3  I'ANDaK!) 
JtW  .An  IN 


V  .  /  )«  * 

0.6727 

0.5413  DEPENDENT  VARIABLE 


NDEPENDENT 
VART  ABLE 


AREA 
L'S**  .  5 


REGRESSION 

COEFFICIENT 


0.780995 

-0.034300 


PARTIAL 

DETERMINATION 

COEFFICIENT 

0.9004 

0.0000 


REGRESSION 
CONSTANT 
1 . 391 2S1 


R 

SQUARE 

0.9993 


R  BAR 
SQUARE 
0.9930 


STANDARD 
ERROR  OF 
ESTIMATE 
0.0705 


VARIABLE  DELETED  IS  L/S**.5 


INDEPENDENT 

VARIABLE 

REGRESSION 

COEFFICIENT 

PARTIAL 

DETERMINATION 

COEFFICIENT 

AREA 

0.749774 

0.9972 

REGRESSION 

CONSTANT 

1.413246 

R 

SQUARE 

0.9997 

STANDARD 

R  BAR  ERROR  OF 

SQUARE  ESTIMATE 

0.9972  0.0614 
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******************************** 

MULTIPLE  LINEAR  REGRESSION 

704-G 1-L2020  JAN  1975 

REVISED  JAN  1933 

******************************** 


TEST  NO.  3  JANUARY  1975 
USE  DATA  FROM  TEST  NO.  1 
READ  REGRESSION  COEFFICIENTS 


NVAL 

NCOMB 

IPRNT 

IFRMT 

I  DELE 

DELTA 

4 

0 

1 

0 

0 

0.0000 

+  +  *■  1-  ANALYSIS  NO  1  +  P  f  +• 
DEPENDENT  VARIABLE  —  LOG  0 


IDEP  . TRANSFORMATION  CODES .  NOBR  IRES  IFORC 

1  1  1  1  100  0  00000000000  0  -1  0 

VARIABLE  TRANSFORMATION 
LOG  Q  NONE 

LOG  3 NO  NONE 

LOG  GW  NONE 

LOG  PRCP  NONE 


STATISTICS  OF  DATA 


BASED  ON  12 

OBSERVATIONS 

STANDARD 

VARIABLE 

AVERAGE 

VARIANCE 

DEVIATION 

LOG  SNO 

0 . 2960 

0.0050 

0.0704 

LOG  GW 

0.4005 

0.002B 

0.0531 

LOG  PRCP 

0.9421 

0.0572 

0.2392 

LOG  Q 

0.9196 

0.0181 

0.1346 

DEPENDENT  VARIABLE 

STMPLF  CORRELATION  COEFFICIENTS 


VARIABLE 

LOG  SNO 

LOG  GW 

LOG  PRCP 

LOG  O 

LOG  SNO 

1.0000 

0.0000 

-0 . 0459 

0.6309 

LOG  GW 

0.0000 

1  .  0000 

0. 1275 

0.4170 

LOG  PRCP 

-0 . 0459 

0.1275 

1 . 0000 

0.2011 

LOG  g 

0  C'3  09 

0  4170 

0  201  1 

1  0000 

IN  ol  11 


INPUT  REGRESSION  COEFFICIENTS 


CONSTANT 


-0.2200 


LOG  3NO  1.6000 

LOG  GW  1.0000 

LOG  PRCP  0.3000 


DEPENDENT  VARIABLE  LOG  Q 


ERROR  AS 


OBS  NO 

OBS  ID 

OBSERVED 

CALCULATED 

RESIDUAL 

DEVIATE 

RATIO 

1 

1  336 

0.939 

0.956 

-0.017 

-0.129 

-0.019 

2 

193  7 

0 . 94  5 

0.904 

0.04  1 

0.  305 

0.043 

3 

1939 

1.052 

1 . 044 

9.008 

0.058 

0.009 

4 

1939 

0.744 

0.776 

-0.032 

-0.237 

-0.043 

5 

194  0 

0.666 

0 . 694 

-0.028 

-0.206 

-0.042 

6 

1941 

1 . 091 

1.110 

-0.029 

-0.213 

-0.027 

7 

194  2 

1 . 060 

1.098 

-0.038 

-0.286 

-0.036 

9 

194  3 

0 . 992 

0.937 

-0.045 

-0.338 

-0.051 

9 

1944 

1.021 

1.019 

0.002 

0.015 

0.002 

10 

194  5 

0.920 

0 . 997 

-0.077 

-0.571 

-0.084 

1 1 

1  346 

0.75  5 

0.773 

-0.018 

-0.135 

-0.024 

1  2 

194  7 

0.960 

0.932 

0.028 

0.210 

0.029 

R  BAR  SQUARED 
STANDARD  ERROR  OF  ESTIMATE 
MEAN  ERROR 
MEAN  SQUARED  ERROR 


0.9190 

0.0383 

-0.0171 

0.0013 
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DESCRIPTION  OF  INPUT  DATA 
GENERALIZED  COMPUTER  PROGRAM 
704-G1 -L2020 


T1 ,  T2 ,  T3  Title  Cards 

Three  title  cards  for  the  identification  of  the  ouput.  A  T1  must  be 
in  columns  1  and  2  of  the  first  card.  A  new  T3  card  is  required  for 
each  subsequent  analysis  which  uses  the  same  data  set. 


J1  -  First  Job  Card 

Field^  Variable  Value 


NVAL 


2.  NCOMB 


3.  IPRNT 

0 

1 


Description 

Total  number  of  values  (variables)  to 
be  read  for  each  observation  and  stored 
for  use  in  this  and  subsequent  analyses. 
This  is  not  necessarily  the  number  of 
variables  to  be  included  in  the  first 
analysis  which  will  be  specified  by  the 
TR  card.  Dimensioned  for  18  variables. 

Number  of  new  variables  to  be  computed 
by  combining  two  or  more  variables  and 
constants.  All  combinations  for  subse¬ 
quent  analyses  which  use  the  same  data 
must  be  performed  in  first  analysis. 
Dimensioned  for  five  combinations.  Sum 
of  NVAL  and  NC^MB  must  not  exceed  18. 

Data  list  options . 

Provides  list  of  input  data  for  each  job. 
Suppresses  listing  of  input  data. 


(1) 


I 


The  standard  card  format  for  most  HEC  programs  consists  of  columns  1 
2  reserved  for  the  card  code  (Field  0)  and  the  remaining  columns  used 
for  program  specification  or  data.  Columns  3-8  are  Field  1  and  the 
remaining  72  columns  are  divided  into  nine  8-column  fields. 
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J1  Card  (continued) 


Field  Variable  Value 
4 .  I FRMT 


0 


1 


5.  IDELE 


0 


1 


6-10. 


NM_  -  Variable  Name  Card 

Field  Varjab'  e  Value 

1 .  DELTA 


2-10.  ANAMA 


Description 

Format  of  data  on  the  cards. 

Data  in  standard  format  of  (A2 ,A6 ,9F8 . 0) 
for  card  identification,  observation 
identification,  and  data,  respectively. 

The  format  for  the  second  card,  if 
required,  is  (8X,9F8.0). 

Special  format  will  be  read  on  FT  card. 

Deletion  options. 

Automatic  deletion  of  each  variable 
with  lowest  partial  determination  coefficient, 
in  turn,  until  only  one  independent  variable 
remains . 

Suppresses  automatic  deletion. 

Not  used. 


Description 

Increment  which  will  be  added  to  al 1  data 
to  meet  certain  constraints  when  logarithmic, 
square  root,  or  reciprocal  transformations 
are  made.  The  increment  should  be  zero 
unless  the  following  constraints  are 
violated: 

Transform  Constrai nt 

Logarithmic  "data"  +  DElTA  >0 

Square  root  "data"  +  DELTA  >0 

Reciprocal  "data"  +  DELTA  -9 

Alphanumeric  name  of  each  variable  read 
in,  NVAL  (Jl.l)  names.  Order  must  corres¬ 
pond  to  order  of  data  as  read.  If  more 
than  one  card  is  needed,  the  name  for  the 
10th  variable  is  placed  in  field  2  of 
the  second  card,  etc. 
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FT  -  Special  Format  Card 

Supply  thib  card  only  if  IFRMT  (J1.4)  is  positive. 

Variable  Val ue  Description 

IFMT  Format  of  data  to  be  read  from  cards. 

The  first  field  must  read  the  data  card 
code  (the  DT)  in_AT_format  and  the  second 
field  must  provide  for  an  alphanumeric 
identification  of  the  observation  in 
A  format  and  NVAL  (Jl.l)  floating  point 
fields,  e.g.,  (A2,  A4,  10X,  10F6.0).  The 
second  A  field  cannot  exceed  six  columns. 
NOTE :  Parentheses  must  be  provided  in 
specification. 

/  t  \ 

CJ  -  Combination  Job  Card1,1 


Field 

1-10. 


Field  Variable  Value  Description 


1 . 

NOP 

Number  of  operation  codes  necessary  to 
compute  a  new  variable  by  combining  one 
or  more  variables.  Dimensioned  for  20 
operation  codes. 

2. 

NC0N 

Number  of  constants  necessary  to  compute 
a  new  variable.  Dimensioned  for  five 

constants . 

3. 

ANAMA 

Alphanumeric  name  of  new  variable  being 
computed . 

4-10. 

Not  used. 

C 0  -  Combination  Operations  Card ^ ^  ^ 

Supply  this  card  only  in  conjunction  with  CJ  card. 

Field  Variable  Value  Description 

1-10-  I0P  Operation  to  be  performed,  NOP  (CJ.l) 

values. 


f 


TTT 


Provide  NCOMB  (J1.2)  sets  of  CJ,  CO,  and,  if  required,  CC  cards. 
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CO  Card  (continued) 


The  sequence  of  operations  necessary  to  obtain  a  desired  computed 
variable  must  be  specified  by  "operation  codes."  These  operations 
are  performed  by  using  three  operating  registers,  argument  (ARG), 
accumulator  (ACC)  and  storage  (STORE).  The  accumulator  has  been 
initialized  to  zero  prior  to  computations .  Any  constants  necessary 
for  the  arithmetic  will  be  entered  on  the  CC  card  in  the  sequence 
necessary.  The  final  value  in  the  accumulator  is  used  as  the  value 
for  the  new  variable. 


Operation 

Var  iable->ARG  -vn 

ARG+ACC-ACC  1 

ACC-ARG-»ACC  2 

ARG  (ACC) -“-ACC  3 

ACC/ARG-ACC  4 

ACCARG-ACC  5 

ACC+STCRE 

then  0.->-ACC  6 

STCRE-ARG  7 

Corstant-^ARG  8 


CC  -  Combination  Constants  Card^  ^ 


Operation  Code 

Variable  number,  preceded  by  minus 
sign,  to  use  next  where  "vn"  is  order 
number  as  read  in  or  combined.  This 
variable  becomes  the  argument. 

Add  argument  to  accumulator. 

Subtract  argument  from  accumulator. 

Multiply  argument  times  accumulator. 

Divide  accumulator  by  argument. 

Raise  accumulator  to  power  of  argument. 


Store  accumulator  while  further  arith¬ 
metic  is  done.  Accumulator  is  then 
reset  to  zero. 

Retrieve  stored  value.  This  value  now 
becomes  the  argument. 

Use  the  next  constant  on  the  CC  card 
as  the  argument. 


Supply  this  card  only  if  NCON  (CJ.2)  is  positive. 

Field  Variable  Val ue  Description 

^ _5  _  Cnsj  Supply  NCON  (CJ.2)  values. 

r  in  Not  used. 


^Provide  NCOMB  (J1 .2)  sets  of  CJ ,  CO,  and,  if  required, 


CC  cards. 
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-  Analysis  Specification  Card 

This  card  is  provided  for  each  analysis. 


Field 


Variable  Value 


Description 


IFUNC 


Number  identifying  the  order  number 
of  the  dependent  variable  in  the  total 
array  of  values  read  in,  i.e.,  3  if 
third  value  in  the  array  is  the  desired 
dependent  variable. 

Variable  deletion  or  transformation  code 
for  each  variable.  NVAL  (Jl.l)  plus 
NCOMB  ( J 1 . 2 )  values  in  special  (12) 
format,  i.e.,  the  code  for  the  first 
variable  read  would  be  in  column  10,  the 
code  for  the  second  variable  in  column  12, 
etc. 


Operation 

Delete  from  current  analysis 
No  transformation 


Square  root 
Logarithmic,  base  10 
Reci pr jcal 

The  number  of  observations  to  be  used 
in  this  analysis.  The  program  will 
count  the  number  of  observations  (NOBS) 
as  the  data  are  read;  therefore,  may  be 
left  blank  if  all  observations  are  to 
be  used  in  analysis.  If  positive,  the 
first  NOBR  observations  will  be  used, 
but  must  be  less  than  or  equal  to  NOBS. 
Dimensioned  for  500  observations. 

Residual -prediction  options.  Provides 
a  tabulation  of  the  observed,  calculated, 
the  difference  between  the  observed  and 
calculated,  and,  if  the  log  transformation 
of  the  dependent  variable,  the  ratio  of 
the  observed  to  the  calculated.  The 
computation  for  a  multiple  linear  regres¬ 
sion  analysis  will  be  performed  when  IRES>0. 
Predetermined  regression  parameters  will 
be  read  (RP  card)  when  I RES=  -1. 
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TR _ C ard  (continued) 

Field  Variabl e  Value  Description 

-1  Regression  parameters  to  be  read  in 

with  no  analyses . 

0  Suppresses  residual -- prediction 

routine . 

>0  Positive  number  indicating  maximum 

number  of  independent  variables  to 
be  included  in  equation.  Routine 
will  operate  for  each  deletion  there¬ 
after. 

9-10-  Not  used. 

RP  -  Regression  Parameters  Card 


Supplv  this  card  only  if  IRES  (TR.8)  is  -1 . 


Field 

Variable 

Val  ue 

Description 

1 . 

AA 

Regression  constant. 

2. 

BB 

Regression  coefficients.  Supply  same 
number  of  coefficients  as  there  are 
independent  variables  (one  less  than 
the  number  of  IFUNC  (TR.2)  values 
that  are  greater  than  zero)  in  the 
order  corresponding  to  data. 

DT  -  Input  Data  Cards 


Input  data  need  only  be  provided  for  first  analysis.  Subsequent  analyses 
can  use  same  data  provided  all  necessary  data  are  read  in  first  analysis. 
The  data  may  be  in  the  following  format  or  the  format  specified  by  the  FT 


card. 

Dimensioned  for 

500  observations . 

Field 

Variable 

Val  ue 

Description 

1 . 

XID 

An  alphanumeric  identification  for  the 
observation. 

2. 

X 

Value  for  the  first  variable. 
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DT  Card  (continued) 


Field  Variable  Value  Description 

3_10.  X  Value  for  succeeding  variables,  NVAL 

(Jl.l)  values.  If  more  than  one 
card  is  needed,  the  value  for  the 
10th  variable  is  placed  in  field  2 
of  the  second  card,  etc.  The  identifi¬ 
cation  (XID)  can  be  repeated  in  field  1 
to  assist  in  identifying  each  card. 

ED  -  End  of  Data  Card 


Supply  tne  same  number  of  ED  cards  as  required  to  supply  one  observation 
on  UT  cards. 


Field  Variable  Value 

0.  ID 

1-10. 

NJ  -  Next  Job  Card 

Field  Variable  Value 

1 .  NEXT 


Description 

The  letters  ED  must  appear  in  columns  1 
and  2. 

Not  used. 


Description 

This  variable  indicates  what  type  of 
job  follows  according  to  the  following 
values : 


0  Begin  a  new  job.  Branches  to  beginning 

of  program  and  reads  TI ,  T2,  and  T3 
cards  for  a  new  job  or  for  a  normal 
stop. 


1  Different  analysis  with  same  input 

data.  Branches  to  read  T3,  TR,  and, 
if  required,  RP  cards.  (No  data  cards 
provided . ) 


Tl  card,  plus  three  blank  cards,  will  cause  a  normal  STOP. 
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/  Either  T3  or  T1  depending  on  NEXT  (NJ.l).  T1  card  plus 

(. _ J  _ A  -thr^e  blank .cards  M.i.l 1  cause_-a_  normaJ  SIQE 

NJ^/  NEXT  |  |  |  |  |  |  ! 

ED/  'l 

DT  /  XID  X  X  ...  NVAL  val ues  .Standard  format  unless 

( _ 1  1 _ I  I  1-1 _ |  IFRMT  ( J1  .4) >0 

RP  j/  AA  BB  BB  ...  NINDP  values.  Omit  unless  IRES  ) 

/  _J _ l _ __1 _ 1 _ 1 _ _L _ J _ 1TR,8JL>0 _ 

TrKTDEP  IFUNC  ..  NVAL  *  NCPf'T  values  NOBR  IRES  ) 

/ _ | _  Special  12  format  i _ L___J _ J _ 

CC  / CNST  CNST  ....  NCON  values.  Omit  unless  NCON  (CJ.2)>0  / 

L _ I _ J _ i _ J _ 1 _ L _ i___J _ 1 _ _ 

( _ _1 _ 1  __  _1_J _ L _ I. _ J _ L _ J _ 

_ L _ l _ J _ 1 _ L _ jLul.2)>h _ ^ 

FT /(Nonstandard  FORMAT  of  input  data).  Omit  unless  IFRMT  |  , 

L _ l  J _ I _ I _ L  1 _ ]{JlJjLt£L_ 

£ _ 1 _ I  .  -  1 _ ! _ I _ I _ 1 _ L _ ! _ _ 

J1  /NVAL  NCOMB  IPRNT  IFRMT  IDELE  )  J 

£—3 _ ! _ I _ 1 _ l _ I _ J _ 1 _ J _ N 

T3 /TITLE  CARD  )  J 

T2 /TITLE  CARD  ,  )  ) 

L _ 1 _ i _ j _ j _ i _ L _ i  i  j _ , 

T1  /tYtTe  card  ^ 
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